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^ ■ Abstract 
(N 

in 

, origin was usually proposed as the most feasible mechanism for the acceleration, there appeared some gen- 

■ eralizations of Einstein equations which could mimic dark energy. In this work we investigate observational 



Recent type la supernova measurements and other astronomical observations suggest that our universe 
is, at the present epoch, in an accelerating phase of evolution. While a dark energy of unknown form and 



constraints on a modified Friedmann equation obtained from the generalized Lagrangian C oc i?" minimally 
coupled with matter via the Palatini first-order formalism. We mainly concentrate on such restrictions of 
model parameters which can be derived from distant supernovae and baryon oscillation tests. We obtain 



o 
U 

■ confidence levels for two parameters (n, flm,o) a-nd find, from combined analysis, that the preferred value 

of r^m.o equals 0.3. For deeper statistical analysis and for comparison of our model with predictions of the 
' ACDM concordance model one applies Akaike and Bayesian information criteria of model selection. Fi- 

5^ , nally, we conclude that the Friedmann-Robertson- Walker (FRW) model merged with a first-order non-linear 

gravity survives SNIa and baryon oscillation tests. 
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I. INTRODUCTION 



The recent observations of type la distant supernovae indicate that our Universe is currently 
accelerating P, Q|. There are different proposals for explaining this phenomenon. Some of them 
are based on assumptions of standard cosmological models, which utilize FRW metric. Thus 
)ossible explanations include: cosmological constant A 0,0], a decaying vacuum energy density 
^ , an evolving scalar field or quintessence models , a phantom energy (expressed in terms of the 
barotropic equation of state violating the weak energy condition) PjQ) dark energy in the form 
of Chaplygin gas P], etc.. All these conceptions propose some kind of new matter of unknown 
origin which violate the strong energy condition. The Universe is currently accelerating due to the 
presence of these dark energy components. 

On the other hand, there are alternative ideas of explanation, in which instead of dark energy 
some modifications of Friedmann's equation are proposed at the very beginning. In these ap- 
proaches some effects arising from new physics like brane cosmologies, quantum effects, anisotro^ 
effects etc. can mimic dark energy by a modification of Friedmann equation. Freese & Lewis 
have shown that contributions of type p" to Friedmann's equation 3H^ = Peff, where p is the 
energy density and n a constant, may describe such situations phenomenologically. These models 
(by their authors called the Cardassian models) give rise to acceleration, although the universe 
is flat and contains the usual matter and radiation without any dark energy components. In the 
authors' opinion 10], what is still lacking is a fundamental theory (like general relativity) from 
which these models can be derived after postulating Robertson Walker (R-W) symmetry. We ar- 
gue that a possible candidate for such a fundamental theory can be provided by non-linear gravity 
theories (for a recent review see e.g. and references therein) and, particularly, the so-called 
f{R)- theories [l^ . It is worth pointing out that if one imposes the energy-momentum conservation 
condition then matter density is parametrized by the scale factor (in a case of R-W symmetry), 
and Cardassian term in the Friedmann equation will be reproduced. 

There are different theoretical attempts to modify gravity in order to achieve an accelerating 
cosmic expansion at the present epoch. Already in the paper by Carroll et al. one can find 
interesting modifications of the Einstein-Hilbert action with Lagrangian density £ oc R+f{R, P, Q). 
Those authors have shown that in the generic case cosmological models admit, at late time, a de 
Sitter solution, which is unfortunately unstable. Moreover, Carroll et al. have demonstrated the 
existence of an interesting set of attractors, which seem to be important in the context of the dark 
energy problem. 
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The main goal of the present paper is to set up observational constraints on parameters of 
cosmological models inspired by non-linear gravity. The possibility of explaning cosmic acceleration 
in terms of nonlinear generalization of the Einstein equation has been previously addressed in 
14,ll5||. However, these authors have not confronted their models with observational data. This 

n 

problem has been tackled in [la], where nonlinear power law lagrangians were compared with SNIa 
data and X ray gas mass fraction as well (see also for a more general class of lagrangians). 
Here, we use samples of supernovae la 18,^^ together with the baryon oscillation test for 
stringent and deeper constraint on model parameters. We check to which extent the predictions 
of our model are consistent with the current observational data. 

Severe constraints on the particular modifications of gravity considered in this paper have been 



already proposed 
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le 



22 



r^, [^. On the other hand, in the article by Clifton and Barrow 
strong constraints coming from nucleosynthesis of light elements have been found within a higher- 
order gravity. Therefore, it is possible that our model (although a first-order), which fits SNIa data 
well can be ruled out by nucleosynthesis arguments. 

Assuming FRW dynamics in which dark energy is present, the basic equation determining a 
cosmic evolution has the form of a generalized Friedmann equation 



tt2 _ PcS k /-|\ 



where Pesio.) stands for the effective energy density of several "fluids", parametrized by the scale 
factor a, while k = ±1,0 denotes the spatial curvature index. One can reformulate in terms of 
density parameters as 

—^ = nes{z) + nk,o{i + zf, (2) 

where ^ = i^csiz) = ilm,o(l + z)^ + ^x,of{z) and Om,o is the density parameter for the 
(baryonic and dark) matter, scaling like a~^, while f{z) describes the dark energy X. For a = uq 
(the present value of the scale factor which we further on normalize to unity), one obtains the 
constraint ^eS,o + ^fc,o = 1- 

We can certainly assume that the energy density (i = m, X) satisfies the conservation condition 



Pi = -3H{pi+pi) (3) 

for each component of the fluid, so that p^s = Spj. Then from Q one gets the constraint relation 
SjOj^o + ^fc,o = 1 foi' the present values (z = 0) of the density parameters. 
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All approaches mentioned above lead toward a description of dark energy in the framework 
of standard FRW cosmology. It will be demonstrated, in the next Section, that all cosmological 
models of the first-order non-linear gravity which satisfy R-W symmetry, can also be reduced to 
the familiar form ((JJ. Therefore, the effects of nonlinear gravity can mimic dynamical effects of 
dark energy. 



II. FRW COSMOLOGY AND FIRST-ORDER NON-LINEAR GRAVITY 



For the cosmological applications one chooses the Priedmann-Robertson- Walker metric, which (in 
spherical coordinates) takes the standard form: 

1 



g = -dr + a\t) 



(4) 



1 — /^r^ 

As before, a{t) denotes the scale factor and k the spatial curvature {k = 0, 1, —1). Another main 
ingredient of all cosmological models is a perfect fluid stress-energy tensor, expressed by 

/ p n n \ 







1—kr'^ 












(5) 





pa^{t)r^ 
pa'^{t)r'^sm'^{d) J 

One requires the standard relations between the pressure p, the matter density p, the equation of 
state parameter w and the expansion factor a{t), namely 



p = wp , p = rja 



Let us consider the action functional 



-3{l+«>) 



7] = const. 



A = ^grav + ^mat = j (x/dit^/(i?) + 2KL^^t{^))d' 



X 



(6) 



(7) 



within the first order Palatini formalism 15]. In fact, from now on we shall assume the simplest 
power law gravitational Lagrangian of the form 



n 



where one fixes the constant 13 to be positive (it has the same dimension as E} We want to 

n 

point out that our model is singular for n = 2. As shown in [1^, such models are exactly solvable 
for the matter stress-energy tensor representing a single perfect fluid of a kind w (cf. (El)). Their 



confrontation with experimental data has been performed, for a dust filled universe, in 



lq |. Here 



we attempt to continue the analysis with newly available Astier SNIa samples and new baryon 
oscillation tests. These allow us to strengthen the admissible constraints on model parameters. 
Moreover, we extend our research to a matter stress-energy tensor containing two components, 



both with p = wp: a perfect fluid w = const 7^ i and a radiation 



w 



^. It turns out that the 



presence of the radiation term crucially changes the dynamics of our model at the early stage of 
its evolution. In addition, as to be demonstrated in Section III, although one cannot obtain better 
constraints from SNIa data, the combined analysis of SNIa and baryon oscillations offers a new 
possibility for a deeper determination of model parameters. 

n 

Following a method developed in the Hubble parameter for our model can be calculated 
to be: 



2n 



3{3w - l)[3w{n - 1) + (n - 3)] 



k{1 — 'iw)rj^ 



/3 



+ 



4n(2 - n)Kr]rad 



2,(3\2,w{n - 1) + (n - 3)]^ 



k{1 — 3w)riu 



n + 3 + 3 i(J ( 1 — ri ) 



k 



2n 



3(l+m) 

a " + 
n 2 



3w{n — 1) + (n — 3) 



(8) 



(Since radiation is already included in (jHJ, one has to assume w 7^ |.) 

It is worth pointing out that the deceleration parameter, in the case k = "q^ad = 0, equals to 
(see Q): 



q{n,w) 



3(1 + - 2n 
2^ 



-1 + 



3(1 + w) 
2n 



(9) 



Thus, the effective equation of state parameter Weff is 

1 + w 



eff 



-1 + 



n 



(10) 



Let us observe that, in the case ry^ad 7^ 0, w < |, the same values of q{n,w) and Weff can also 
be achieved as asymptotic values a 00. In the early universe, when the scale factor goes to the 
initial singularity, the radiation term, in (jH)) scaling like a~^^^"\ will dominate over the dust term 
(scaling like a^~). More precisely, if n < or n > 2, then the negative radiation term cannot 
dominate over the matter, so that instead of the initial singularity we obtain a bounce. On the 
other hand, if a goes to infinity, the radiation becomes negligible versus to the matter. 
In our further analysis we restrict ourselves to the case w = k = i.e., more precisely, to the 
spatially flat universe filled with dust and radiation. Thus (8), (remarking once more that all r/'s 
are positive constants) becomes 



2n 
3(3 - n) 



_ 3 4n(2 - n)Kr]rad 
a " + 



3/3(n - 3)2 



Vdustf^ 



n+3 

a " 



(n - 3)2 



(11) 
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One should immediately note that this expression, representing the squared Hubble parameter, 
reproduces in the case n = /3 = 1 - as expected - the standard Priedmann equation. We would 
like to emphasize also that (11) becomes singular at n = 3. 



It is convenient to rewrite relation (11) is such a way that all coefficients are dimensionless 
(density parameters). Then, the effects of the matter scaling like a"^^^"*"""^, and the radiation 
scaling like a~^, can be separated from the effects of the nonlinear generalization of Einstein 
gravity (n / 1): 



H\ N'^ 2n ^ ,^ , 3(l-n) 

Ho J (3-n) 

+0,,o(l + Z)^ ^'^^'^~^ ^nonl,o{^ + z)^^ . (12) 

(n-3) 

Here: flm,o = ^^f^T^, ^r,o = ^^^i ; ^nonifi = ^^d^t^^ " while Hq denotes the present-day value 
of the Hubble function. Let us observe that i^nonifl can be determined also from the constraint 
H{z = 0) = Hq, which easily reduces to: 

2n _ 4n(2-n) ^ 



V(3-n) (n-3)^ 
The relation (|12|) has the form of Friedmann's first integral. Therefore, the dynamics of the 
model can be naturally rewritten in terms of a 2D dynamical system of Newtonian type. Its 
Hamiltonian is: 



7^ = 1^2 ^ y{a) = , (13) 
while the corresponding equations of motion are: 



a = 2/ , 

The overdot differentiats now with respect to the rescaled time variable r, so that dt = \KQ\dT, 
while V{a) is a potential function for the scale factor a expressed in units of its present value 
ao = 1. If = /(a), then the potential function is given by the general formula 



^(a) = -\f{aV ■ (15) 



For example, the potential function for our model writes as (ag = 1): 



2n _i 4n (2 - n) _2 

\ "m.OQ^ H -Q-^KfiO, 

> — n- (^72 — 3) 



f^noni.oa-^'-"^ ■ (16) 



Recently, in Carloni et al .2^, the cosmological dynamics of i?" gravity has been treated in a 
different phase space with the use of qualitative methods for dynamical systems. 

The phase portraits for the ACDM model versus our model with fitted values of n, i^mfl pa- 
rameters (see the next Section) are illustrated on Figures 1 and 2. Both models are topologically 
inequivalent: the phase portrait of ACDM has a structurally stable saddle critical point, while 
with nonlinear gravity one obtains a center. As well-known, the critical point of a center type is 
structurally unstable and all trajectories around this point represent the models, which oscillate 
without initial and final singularities. 

——. 3(l-n) 

It is interesting that H12|) after a time reparametrization following the rule: dr/ = (1 + z) 2 dr 
is equivalent to the standard cosmological model with matter and radiation, with rescaled values 
of the corresponding density parameters r^m.o and ilrfi. 

The geometry of the potential function offers the possibility to investigate the remaining models. 
On can simply establish some general relation between the geometry of the potential function and 
critical points of the Newtonian systems. In any case, the critical points lie on the a-axis, i.e. they 
represent the static solution yo = 0, a = qq so that (^)^^ = 0. If (oq, 0) is a strict local maximum 
of V{a), it is of the saddle type. If (ao, 0) is a strict local minimum of the analytical function V{a), 
it is a centre. If (oq, 0) is a horizontal inflection point of the V{a), it is a cusp. 

From the fitting procedure we obtain n > 2, so the second term in the potential function is 
negative (in contrast to the first term which is positive). Because the negative radiation term in 
(|16|) can not dominate the first one {V <0), there is the characteristic bounce behavior rather than 
the initial singularity in the ACDM model. Moreover, during the bouncing phase the universe is 
accelerating, while for late times it becomes matter dominated and decelerates. 



III. DISTANT SUPERNOVAE AS COSMOLOGICAL TEST 

Type la distant supernova surveys suggest that the present Universe is accelerating P, O] . Every 
year new SNIa enlarge the available data by more distant objects and lower systematics errors. 
Riess et al. 3] have compiled samples which become the standard data sets of SNIa. One of 
them, the restricted "Gold" sample of 157 SNIa, is used in our analysis. Recently Astier et al. [l^ 
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FIG. 1: The phase portrait for the ACDM model. There is a single critical saddle-point on the a-axis. It 
represents the static Einstein universe. The trajectory of flat fc = model divides all remaining ones into 
closed (inside) and open (outside) models. 

have compiled a new sample of supernovae, based on 71 high redshifted SNIa discovered during 
the first year of Supernovae Legacy Survey. This latest sample of 115 supernovae is used as our 
basic sample. 

For distant SNIa one can directly observe their apparent magnitude m and redshift z. Because 
the absolute magnitude M is related to the absolute luminosity L, the relation between luminosity 
distance di, the observed magnitude m and the absolute magnitude M has the following form 

m - M = SlogioC^L + 25 . (17) 
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FIG. 2: The phase portrait for nonhnear gravity with £ cx R", n = 2.6 (from estimation). There is a single 
critical point on the a-axis - a center; a — a^rit — ^^(^^^^yi The trajectories of the system lie in 

the physical region {a : a > Ssanj ^nd represent bouncing evolution. In this scenario, instead of the big 
bang singularity of the ACDM model, one has a bounce a — amin, a > 0. It lies in a neighborhood of the 
minimum of the potential function. During the bounce phase, the universe is still accelerating. Note that 
if radiation effects vanish, there is no static critical point on the a-axis (formally acrit = is allowed for 
n > 3). 

It is convenient to use the dimensionless parameter Dl 

Dl = HodL (18) 
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instead of di. Then (fT7|l can be replaced by 

lj, = m- M = 5 logio Dl + M , (19) 

where 

7W = -51ogio^o + 25 . (20) 

We know the absolute magnitude of SNIa from its light curve. The luminosity distance of super- 
novae is a given function of the redshift: 

^ ^ '^i^l^ (^o^f . pi) 

where Qk,o = and 

T{x) = sinh(3;) for k < , 

T{x) =x for = , (22) 

J^{x) = sin(x) for k > . 

Substituting (PT|) back into equations ((T7|) and (fTU)) provides us with an effective tool (Hubble 
diagram) to test cosmological models and to constrain their parameters. Assuming that supernovae 
measurements come with uncorrelated Gaussian errors, one can determine the likelihood function 
C from chi-square statistic £ oc exp(— x^/2), where 

/..theor ,,obs'\2 

X^ = E ^^^ ^2^' ^ • (23) 

i * _ 

The Probability Density Function (PDF in short) of cosmological parameters [l| can be derived 
from Bayes' theorem. Therefore, one can estimate model parameters by using a minimization 
procedure. It is based on the likelihood function as well as on the best fit method minimizing x^. 

For statistical analysis we have restricted the parameter r^m.o to the interval [0, 1] and n to 
[—10.0, 10.0] (except n = and additionally n = 3 for w = 0). Moreover, because of the singularity 
at n = 3,tt; = (see eq. (|12() ) we have separated the cases n > 3 and n < 3 for w = in our 
analysis. Please note that O,noni,o is obtained from the constraint H{z = 0) = Hq. 

In Figure 3 we present residual plots of redshift-magnitude relations (Hubble diagram) between 
the Einstein-de Sitter model (represented by zero line) and our best fitted model — upper curve — 
and ACDM model — middle curve. One can observe systematic deviations between these models 
at higher redshifts. The non- linear gravity model predicts that high redshifted supernovae should 
be fainter than those predicted by the ACDM model. 
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FIG. 3: Residuals (in mag) between the Einstein-de Sitter model (zero line), the flat ACDM model (middle 
curve) and the non- linear gravity model (upper curve). Results obtained with the Astier (left panel) and 
the Riess (right panel) samples. 
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FIG. 4: PDF obtained with the Astier sample for the parameters flm,o and n, marginalised over the rest of 
parameters. Non-linear gravity model (u; = 0, n < 3). 

The results of two fitting procedures performed on Riess and Astier samples with different prior 
assumptions for n are presented in Tables 1 and 2. In the Table 1 the values of model parameters 
obtained from the minimum of given, whereas in Table 2 the results from marginalised 

probability density functions are displayed. Please note that we obtained different values of Ai 
from the Riess versus Astier samples. It is because Astier et al. assume the absolute magnitude 
M = —19.31 ± 0.03 + Blogi^hiQ 19]. For comparison we present (Table 3) results of statistical 
analysis for the ACDM concordance model. 

The best fit (minimum x^) gives n ~ 2.6 with the Astier et al. sample versus n ~ 2.1 with 
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TABLE I: The flat non-linear gravity model with w = 0. Results of statistical analysis performed on the 
Astier versus the Gold Riess samples of SNIa obtained from best-fit. We separately analysed the case 
n > 3 and n < 3 . 



sample 







M 


Gold 


0.35 


< 0.01 3.001 


15.975 180.7 


n < 3 


0.89 


0.23 2.13 


15.975 181.5 


n > 3 


0.35 


< 0.01 3.001 


15.975 180.7 


Astier 


0.01 


-1.47 3.11 


15.785 108.7 


n < 3 


0.98 


0.08 2.59 15.785 108.9 


n > 3 


0.01 


-1.47 3.11 


15.785 108.7 



TABLE II: The flat non- linear gravity cosmological model {w = 0). The values of the parameters obtained 
from marginal PDFs calculated on the Astier versus the Gold Riess samples. Because of the singularity at 
n = 3 we separately analyze the cases n > 3 and n < 3 . 



sample 


^m,0 




n 


M 


Gold 


0.01 


0.26 


2.11 


15.955 Q 


n< 3 


1.00 


0.26 


2.11 


15.955^0.03 


n > 3 


0.01 


-0.01 


3.001 


15.955^o^o| 


Astier 


0.01 


0.09 


2.56 




n< 3 


1.00 


0.09 


2.56 


15.785+°;°i^ 


n > 3 


0.01 


-0.01 


3.01 


15.785in^ 



the Gold sample. In Figure 4 we present PDF obtained with the Astier sample for the parameters 
ftm,o and n for non-linear gravity model, (case n < 3 marginalised over the rest of parameters). 
Please note that from Fig. 4 we obtain a very weak dependence of PDF on the matter density 
parameter if only Qrn,o > 0.05. 

In Figure 5, confidence levels on the plane {Qm,o,n), for non-linear gravity model, for the case 
n < 3 marginalized over Ai are presented. 

Recently Eisenstein et al. have analyzed baryon oscillation peaks detected in the Sloan Digital 

TABLE III: Results of statistical analysis of the ACDM flat model performed on the Astier versus the Gold 
Riess samples of SNIa as a minimum x^ best- fit. 



sample 


^m,0 


f^A,o M 


Gold 


0.31 


0.69 15.955 175.9 


Astier 


0.26 


0.74 15.775 107.8 
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Sky Survey (SDSS) Luminosity Red Galaxies 



20]. They found 



E{zi)^ \zi^\nk,o\ \^ ' Jo ^(^), 
so that E{z) = H{z)/Hq and zi = 0.35 yield A = 0.469ib0.017. The quoted uncertainty corresponds 
to one standard deviation, where a Gaussian probability distribution has been assumed. These 
constraints could also be used for fitting cosmological parameters |l9l. l24l|. We obtain from this 
test the values of the model parameters ilmfl = 0.28, O,noni,o = 0.33 and n = 2.53 for a best fit. 
In Figure 6 we show the region allowed by the baryon oscillation test on the plane (17^,0;^) for 
non-linear gravity model (for the case n < 3). In Figure 7 we present combined confidence levels, 
obtained from the analysis |3l of both data sets. We find that the model favours J7m,o — 0.3 and 
n ~ 2.6. 

In modern observational cosmology, one encounters the so-called degeneracy problem: many 
models with dramatically different scenarios (big bang or bounce, big-rip or de Sitter phase) agree 



2^ I can be used. 



27l | . Among 



with the present day observational data. Information criteria for model selection 
in some subclass of dark energy models, in order to overcome this degeneracy |2£ 
these, Akaike (AIC) [2^ and Bayesian (BIG) information criteria [29I are the most popular. From 
these criteria one can determine several essential model parameters, providing the preferred fit to 
the data 



The AIG |28|] is defined by 



AIG = -21n£ + 2d , (25) 



where C is the maximum likelihood and d the number of model parameters. The best model, with 
a parameter set providing the preferred fit to the data, is that which minimizes the AIC. 
The BIG introduced by Schwarz [2^ is defined as 

BIG = -21n£ + dln7V , (26) 

where N is the number of data points used in the fit. While AIG tends to favor models with a 
large number of parameters, the BIG penalizes them more strongly, so the later provides a useful 
approximation to the full evidence in the case of no prior on the set of model parameters [s^l . 

The effectiveness of using these criteria in the current cosmological applications has been recently 
demonstrated by Liddle j^l- Analyzing GMBR WMAP satellite data he found the number 
of essential cosmological parameters to be five. Moreover, he came to important conclusion that 

13 




FIG. 5: The flat non-linear gravity model {w = Q,n < 3). Confidence levels on the {n„i,o,n) plane, 
marginalised over A4, obtained from SNIa Astier sample. 

spatially-flat models are statistically preferred to close models as it was indicated by the CMBR 
WMAP analysis (their best-fit value is Qtot,o = "^1^1,0 = 1-02 it 0.02 at la level). 

In the paper by Parkinson et al. 30], also the usefulness of Bayesian model selection criteria 
in the context of testing for double inflation with WMAP was demonstrated. These criteria were 
also used recently by us to show that models with the big-bang scenario are rather preferred over 
the bouncing scenario 32 1. 



14 




FIG. 6: The flat non-linear gravity model {w — 0,n < 3). Confidence levels on the (r^m.o, n) plane obtained 
from baryon oscillation peaks. 

Please note that both information criteria make no absolute sense and only the relative values 
between different models are physically interesting. For the BIC a difference of 2 is treated as 



a positive evidence (6 as a strong evidence) against the model with larger value of BIC 33, |34 1 . 
Therefore one can order all models, which belong to the ensemble of dark energy models, following 
the AIC and BIC values. If we do not find any positive evidence from information criteria, the 
models are treated as identical, while eventually additional parameters are treated as not significant. 
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FIG. 7: The flat non- linear gravity model (w = 0, n < 3). Common confidence levels on the plane (f2m,0: 
obtained from SNIa Astier sample and baryon oscillations. 

Therefore, the information criteria offer a possibility to introduce a relation of weak ordering among 
considered models. 

For comparizing the ACDM and the non-linear gravity models the results of AIC and BIC 
are presented in Tables IIVI Note that for both samples we obtain with AIC and BIC for the 
ACDM model smaller values than for non-linear gravity. We use a Bayesian framework to compare 
the cosmological models, because they automatically penalize models with more parameters to fit 
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TABLE IV: Results of AIC and BIC performed on the Astier versus the Gold Riess samples of SNIa. 

sample AIC BIC 

ACDM Gold 179.9 186.0 

ACDM Astier 111.8 117.3 

Non-Lin.Grav. Gold 186.6 195.8 
Non-Lin.Grav. Astier 114.7 122.9 

the data. Based on these simple information criteria, we find that the SNIa data still favour the 
ACDM model, because under a similar quality of the fit for both models, the ACDM contains fewer 
parameters. 

It is interesting that both models give different predictions for the brightness of the distant 
supernovae (see Fig. 3). The model of modified non-linear gravity predicts that very high redshift 
supernovae should be fainter than predicted by ACDM. So, we can expect future SNIa data to 
allow us discriminating finally between these two models. 

IV. CONCLUSION 

The main subject of our paper has been to confront the simplest class of non-linear gravity 
models versus the observation of distant type la supernovae and the recent detection of the baryon 
acoustic peak in the Sloan Digital Sky Survey data. We find strong constraints on two independent 
model parameters (Qrnfl, n). If we assume n=l, then we obtain the standard Einstein de Sitter 
model filled by both matter and radiation. We estimate model parameters using standard min- 
imization procedure based on the likelihood function as well as the best fit method. For deeper 
statistical analysis, we have used AIC and BIC information criteria of model comparison and selec- 
tion. Our general conclusion is that non-linear gravity fits well (both SNIa and baryon oscillation 
data). In particular we conclude: 

1. Analysis of SNIa Astier data shows that values of statistic are comparable for both ACDM 
and best fitted non-linear gravity model. 

2. The non-linear gravity models with n < 2 can be excluded by combined analysis of both 
SNIa data and baryon oscillation peak detected in the SDSS Luminous Red Galaxy survey 
of Eisenstein at al. |2Q|] on 2a confidence level. 

3. From SNIa data we obtain a weak dependence of the quality of fits on the value of density 
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parameter for matter (Om,o)- However, the combined analysis allowed only value of 0,^,0 
well tuned to its canonical value 0^,0 = 0.3. This value of course is in good agreement with 
present extra-galactic data [s^. 

4. We use the Akaike and Bayesian information criteria for comparison and discrimination 
between the analyzed models. We find these criteria still to favour the ACDM model over 
non-linear gravity, because (under the similar quality of the fit for both models) the ACDM 
model contains one parameter less. 

5. The Hubble diagram implies that very high redshifted supernovae {z > 1.5) should be fainter 
in non-linear gravity model than those predicted by ACDM. So, future SNIa data can allow 
us finally to discriminate between these two models. 

6. The standard general relativity models with n = 1 (without cosmological constant) can be 
excluded by SNIa data on 17a level (as the E-deS model). 

7. The non-linear cosmology can therefore be treated as a serious alternative versus cosmology 
with dark energy of unknown nature. 
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